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Abstract 



Using a nonperturbative method based on asymptotic behaviour 
of Wilson loops we calculate masses of glueballs and corresponding 
I Regge-trajectories. The only input is string tension fixed by meson 

Regge slope, while perturbative contributions to spin splittings are 
defined by standard values. The masses of lowest glueball states 
are in a perfect agreement with lattice results. The leading glueball 
trajectory which is associated with Pomeron is discussed in details 
■ and its mixing with / and /' trajectories is taken into account. 



^ : 1 Introduction 



The problem of existence of glueballs is one of the most interesting in QCD. 
Lattice calculations give definite predictions for a spectrum of such states 
[1-4], but experimental evidences are not conclusive [5,6]. Mixing between 
gluons and gg-pairs makes a separation of glueballs a difficult problem. A 
theoretical study of glueballs in QCD started in [7-10] and is closely related to 
the problem of Pomeron - leading Regge pole, which determines the asymp- 
totic behaviour of scattering amplitudes at very high energies. It is usually 
assumed that Pomeron in QCD is mostly gluonic object [11] and glueball 
resonances with vacuum quantum numbers and largest spins belong to this 
trajectory. Another interesting hypothetical Regge singularity is the "odd- 
eron", which has negative signature and C-parity and can be built out of 
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at least 3 gluons. Most studies of the Pomeron and odderon singularities in 
QCD are based on applications of the perturbation theory [12]. 

In this paper we will address both the problem of spectra of glueballs 
and of the Pomeron (odderon) singularity using the method of Wilson-loop 
path integrals developed in papers [13-15]. The method is based on the 
assumption of the area law for Wilson loops at large distances in QCD, which 
is equivalent to the condition of confinement of quarks and gluons. It has 
been first applied to calculation of spectra of gg-states in [16], baryons in [17] 
and glueballs in [18]. In this latter calculation rotation of the string between 
quarks or gluons was not taken into account, which lead to some distortion 
of mass spectra, notably the Regge slope was l/Scr instead of string slope 
l/27r(T. In this paper we will use for glueballs more accurate calculational 
method developed in [14,15] for gg-state, which yields the correct Regge 
slope (see [9] for numerical data and discussion). We also make a detailed 
study of possible corrections to large distance string dynamics due to small 
distance perturbative gluon exchanges (PGE) and will demonstrate that their 
influence on the mass spectrum of glueball states is rather small and can be 
computed as a correction. This is in contrast to the glueball spectrum in [8] 
where PGE in the form of the adjoint Coulomb potential was assumed as in 
many other papers on the subject. Instead we argue below in the paper, that 
PGE sums up to another series, the BFKL ladder [12], where loop corrections 
strongly suppress the final result , so that PGE can be disregarded in the first 
approximation. Our predictions for masses of lowest spin-averaged glueball 
states in units of are in a perfect agreement with results of recent lattice 
calculations [1-3]. In addition spin-orbit and spin-spin interactions are also 
calculated and found in a good agreement with lattice data. 

The leading glueball Regge-trajectory is calculated in the positive t re- 
gion and is extrapolated to the scattering region of t < 0. The importance of 
mixing among this trajectory and gg-trajectories (/, /') is emphasized; cal- 
culation of these mixing effects yields the leading Pomeron trajectory with 
ap{0) > 1 {ap{0) — 1.1 -=-1.2) in accord with experimental observations [20]. 
An interesting pattern of 3 colliding vacuum trajectories in the region t > is 
observed, which can be important for decay properties of resonances situated 
at these trajectories. 

The paper is organized as follows. The field operators creating glueball 
states in general case and in the background field method are introduced 
in Section 2 and presented in Appendix 1. The general formalism of the 
Wilson loop path integrals and the resulting relativistic Hamiltonian is dis- 
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cussed in Section 2 and Appendix 2. The spectrum of spin-averaged glueball 
states following from this Hamiltonian is obtained in Section 2 and com- 
pared with corresponding lattice results. Spin splittings of glueball masses 
from both nonperturbative and perturbative parts are considered in Section 
3 and Appendix 3 and 4. Resulting glueball spectrum is compared in the Sec- 
tion with lattice calculations. The influence of PGE on glueball masses and 
Regge - trajectories is discussed in Section 4. It is pointed out that effects 
of small-distances on glueball spectrum are small. Three-gluon glueballs are 
considered in Section 5. It is shown that the lowest 3g states have rather 
large mass M^g ^ SAGeV. Relation between glueball Regge trajectories and 
the vacuum Pomeron trajectory is discussed in Section 6. In this Section and 
Appendix 5 effects of mixing between gluonic and gg-Regge-trajectories are 
investigated. 

Possible implications and improvements of the results of the paper are 
discussed in Conclusions. 

2 General Formalism 

Following [13,21], we separate gluonic fields into nonperturbative back- 
ground Bfj^ and perturbative gluons a^, = 5^ + and consider two- 
gluon glueballs, described by the Green's functions 

G,.,Mx,y\x',y') = (vI/(-)(a;,2/)vi/(-*)(x',yO)a,B = 

(r(^")G^^.(x, x')G,,iy, y')r^°"*))B + perm (1) 

where v]>(*")(v]> are glueball operators in the initial (final) state made of 
gluon field and (see Appendix 1 for explicit form of v]/(*")'(°*'*) in lowest 
states); G^^i is the gluon Green's function of field in the background field 
B^, namely 

G^,{x,y) = {x\i-D'5,, - 2igF^,)-'\y)B (2) 

where = — igB^, F^^ is the field strength of the field B^ in the ad- 
joint representation, and averaging over background B^ is implied by angular 
brackets (where subscript B will be omitted). 

-'^Note that the background formaUsm exploiting the 'tHooft identity [21] allows us to 
avoid the double-counting problem and the principle of separation is unimportant provided 
background is characterized by the input string tension a and perturbation is in the 
known as constant 
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Referring the reader for details of derivation to refs. [13,18,21] and Ap- 
pendix 2, one can write for (1) the path integral 



G^^,^,^,{x,y\x\y') = const ds ds'DzDz'e'^-^ {T^'^'^WfT^"'''^) (3) 
where K — j Joi^Ydr, K' is the same with primed z, r, s, and 

(Wf) = trPePFieMw I B^du^ + 2ig f Fdr + 2ig f Fdr'}) (4) 

Here Pb, Pf are ordering operators of the color matrices and F respec- 
tively. As we shall see in section 3, the terms with F generate spin-dependent 
contribution of nonperturbative background, which is calculable and small, 
and we shall treat those terms perturbatively. 

Neglecting F's as a first approximation and omitting for simplicity pro- 
jection operators F^'") , r*^°"*^ which do not influence the form of resulting 
Hamiltonian, one arrives at the Wilson loop in the adjoint representation, 
for which one can use the minimal area law, confirmed by numerous lattice 
data [22] at least up to the distance of the order of 1/m, 

(Wadj) = Z exp{-aadjSmin) (5) 

where we have included in Z self-energy and nonasymptotic corrections, since 
(5) is valid for large loops with size R,T Tg, where Tg is the gluon corre- 
lation length. 

Note, that we could treat (4) by the field correlator method [13,23], keep- 
ing only lowest (Gaussian) correlator {F{x)F{y)). In this case the leading 
term will again have the form (5); we shall use this method to evaluate con- 
tribution of the gluon spin terms F in (4); the result (5) is more general, 
since aadj in (5) contains contribution of all correlators and is not connected 
to the Gaussian approximation. 

Applying now the general method of [14] to the Green's function (3), 
one introduces auxiliary (einbein) function fj,(t) of the real time t instead of 
the proper time s,t, via relation dr — and einbein auxiliarly function 
u{P, t) to get rid of the square root Nambu-Goto form for S^in in (5). As a 
result one defines the hamiltonian H, through the equality G ~ exp{—TH), 
where T is the evolution parameter, taken here to be the center-of-mass time 
T, < t < T. 
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The resulting relativistic Hamiltonian for two spinless gluons looks like 

[14] 

H - I u(f\ I ^^^^^^ I 



+ Z'^^^-' + i/^KA^M/? (6) 



Here /i(t) and z/(/3, t) are positive auxiliary functions which are to be found 
from the extremum condition [14]. Their extremal values are equal to the 
effective gluon energy (/x) and energy density of the adjoint string {v). 

For the case L — Q the extremization over and u yields a simple answer 
[14] , coinciding with the Hamiltonian of the relativistic potential model 

Ho = 2^^ + aadjr (7) 

An approximation made in [16-18] corresponds to the replacement of the 
operators /i(r), i>{t, (3) (which by extremization are expressed through oper- 
ators p, r) by c -numbers, to be found from extremization of eigenvalues of 
Hq. This yields another form, used in [18], 

p2 

Hq = h //o + (^adjr (8) 

A'-o 

as can be seen from Table 3 of ref. [19], eigenvalues of (8) are about 5% 
higher than those of Hq. The value of aadj in (8) can be found from the 
string tension of qq system, since the Casimir scaling found on the lattice 
[22] predicts that 

_ C^jado) _9 

For light quarks the value of a fund is found from the slope of meson Regge 
trajectories and is equal to 

(X/w = ^~0.18GeV^2 (iQ) 

^From that we find 

Gadj ~ 0.40G'eV^' (11) 

In what follows the parameter ^ and its optimal value //q, which enters in 
(8) play very important role. The way they enter spin corrections in Section 3 
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and magnetic moments shows that /Iq plays the role of effective (constituent) 
gluon mass (or constituent quark mass in the equation for the qq system). 

In contrast to the potential models, where the constituent mass of glu- 
ons and quarks is introduced as the fixed input parameter in addition to 
the parameters of the potential, in our approach /iq is calculated from the 
extremum of the eigenvalue of equation (8), which yields 

/^o(n) = V^{^f\Mo{n) = 4/.o(n) 

where a — Uadj for gluons and a — Uf ior massless quarks, and a(n) is the 
eigenvalue of the reduced equation ^ + (a(n) — p — L{L + 1) / p'^)ip = 0. The 
first several values of a{n) and /Xo(n) are given in the Table 1, and will be 
used in Section 3. 

Note that our lowest "constituent gluon mass" ^^{n — L — {)) — Q.b2%GeV 
(for Of — O.lSGey) is not far from the values introduced in the potential 
models, the drastic difference is that po depends on n, L and grows for higher 
states, and is calculable in our case. 

The eigenvalues of i^o, Eq- (7), for L = and different n = 0, 1, 2, ... are 
given in Table 2 for Uadj = \ofund = |0.18Gel^^. The mass spectrum for L > 
is given by eigenvalues of Hq, Eq.(6), and was studied in [14]. With the 5% 
accuracy of WKB approximation one can exploit much simpler expressions 
found in [19], which predict for L > the eigenvalues shown in Table 3. An 
independent numerical estimation of the rotating string spectrum was done 
in [24] and yields similar eigenvalues. 

^Prom Tables 1,2,3 and [24] one can see that mass spectra of the Hamil- 
tonian (6) are described with a good accuracy by a very simple formula 

= L + 2n^ + ci (12) 

27r(7 

where L is the orbital momentum, Ur -radial quantum number and Ci is 
a constant 1.55. It describes an infinite set of linear Regge-trajcctories 
shifted by 2nr from the leading one (n^ = 0). The only difference between 
light quarks and gluons at this stage is the value of cr, which determines the 
mass scale. 

Thus the lowest glueball state with L = 0, = according to Table 2 
and eq. (12) has = 4.04^6^2^ 

It corresponds to a degenerate 0^"*" and 2++ state. 

M = 2mGeV (13) 
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In order to compare our results with the corresponding lattice calculations 
[1-4 ] it is convenient to consider the quantity M / ^/o/, which is not sensitive 
to the choice of string tension a 0. 

/^From these data we have for L = 0, = states the spin averaged mass 

M _ M(0++) + 2M(2++) 1 ^^^^ 



the value 4.61 ±0.1, which should be compared to our prediction M^^^°^[L = 
0,nr = 0)/yo7 = 4.68. 

For radially excited state our theory predicts 

j^theor 

(L = 0,n^ = 1) = 7.0 (15) 



Lattice data [1] give for this quantity 



(L = 0,n^ = 1) = 6.56 ±0.55 (16) 

For L = 1, = 1 states one can define spin-averaged mass in a similar 
way 

M M(0-+) ± 2M(2-+) 1 ^^^^ 



/af 6 ^Of 

lattice data [1-4] yield 



(L = l,n^ = 0) = 6.11 ±0.38 (18) 



which is in a reasonable agreement with our prediction 



-(L = l,n^ = 0) = 6.0. (19) 



For L = 2,nr = and L = 1, = 1 we have the spin averaged state we 
have 

^ ' ' ' 7.0; (L = l,n^ = 1) = 8.0 (20) 



^Note that the value a/ ~ 0.23GeV^ used in lattice calculations differs by about 20% 
from the "experimental" value (10). 
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Lattice data [1] yield respectively 7.7 ± 0.4 and 7.94 ± 0.48. Note that in 
the first multiplet lattice data exist only for 3"'"+. The overall comparison of 
spin-averaged masses computed by us and on the lattice is shown in Table 4. 

Thus we come to the conclusion that the spin-averaged masses obtained 
from purely confining force with relativistic kinematics for valence gluons are 
in a good correspondence with lattice data, which implies that PGE shifts 
of glueball masses in lattice calculations is small. 

3 Spin splittings of glueball masses 

Here we shall treat spin effects in a perturbative way; a glance at our pre- 
dictions in Table 5 and at the lattice results given in Table 6 tells that spin 
splittings in glueball states apart from 2++ — 0"^+ amount to less than 10- 
15% of the total mass, and hence perturbative treatment is justified to this 
accuracy level. 

Before starting actual calculations of spin splittings one should choose 
between two possible strategies (and corresponding physical mechanisms) of 
treating gluon polarizations. In the first approach one insists on the transver- 
sality condition and on the resulting two gluon polarizations as for the free 
gluon [7], (this is the procedure accepted e.g. in [28]). 

In the second approach it is assumed that gluon has acquired nonzero 
mass due to the adjacent string, similarly to the case of W^, where mass 
is created by the Higgs condensate. In this case one has 3 massive gluon 
polarizations and the spin coupling scheme of two gluons can be taken as 
the LS scheme with the characteristic J^*" pattern of lowest levels, which is 
observed in lattice calculations [1-4]. 

Therefore we choose the second approach and consider gluon spin oper- 
ator S^'\i = 1,2, total spin operator S = S^^^ -|- S*^^^ and orbital angular 
momentum L, the total angular momentum J = L -|- S, and assign to each 
level (mass) not only conserved values of J^*", but also values of L, S (which 
in some cases may have admixture of L' = L ± 2, S' = S ± 2, but this 
admixture in general case is small). 

The detailed discussion of the gluon mass generation in the context of 
gauge invariance and symmetry breaking (as also in the electroweak case) is 
relegated to a separate publication. 

The two-gluon mass operator can be written as 

M = Mo(n, L) + SLMsL + S^^^S^^^ Mgs + Mt, (21) 
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where Mq is the eigenvalue of the Hamiltonian H = Hq + AHp^rt, and Hq is 
given in (7) (or its approximation in (8)), while AHpf,rt is due to perturbative 
gluon exchanges and discussed in the next section. 

To obtain three other terms in (21) one should consider averaging of the 
operators F in the exponent of (4) and take into account that 

- 2iF^, = 2(S(^)b(^) + S^'^E(^))^. (22) 

and similarly for the term in the integral / Fdr', with the replacement of 
indices 1 — > 2. Here gluon spin operators are introduced, e.g. 

i^rn^)ik ~ ~'^^miki i,k = 1, 2, 3, (>S'^-')j4 = —iSim (^3) 

Two remarks are in order here: i) the gluon spin enters via the integral 
/ 2SB(ir' = / ^^dt, where fi(t) with its cxtremum value fiQ is the same as in 
(6), (8) (for details see Appendix 2). ii) the main part of the Hamiltonian, Hq 
is diagonal in spin indices i, k; while the spin-dependent part (22) is treated 
as a perturbation, hence the admixture of the 4-th polarization due to S in 
(23) does not appear to the lowest order. 

The detailed derivation of spin-dependent terms is done in the Appendix 
3, here we only quote the results. Since the structure of the term F in (4) 
due to (22) is the same as in case of heavy quark with the replacement of the 
heavy quark mass by the effective gluon parameter /iq (see (8)), one can use 
the spin analysis of a heavy quarkonia done in [25], to represent the spin- 
dependent part of the Hamiltonian in the form similar to that of Eichten and 
Feinberg [26] 



AH - —(-— - — ) ^^^^^V(r) 
* r dr r dr 3/Xq ^ 



+ -^(3(S«n)(S(2)n) - S^^^S^^^)V^{r) + AV (24) 

where AV contains higher cumulant contributions which can be estimated to 
be of the order of 10% of the main term in (24) and will be neglected in what 
follows. Note that spin of gluon is twice that of quark, therefore spin-orbit 
and spin-spin terms for glueballs are effectively twice and four times larger 
respectively than for the quarkonia case. 

The functions Vi{r) are the same as for heavy quarkonia [25] except that 
Casimir operators make them 9/4 times larger; the corresponding expressions 
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of Vi{r) in terms of correlators D{x), Di{x) [23], are given in the Appendix 3. 
Both D and Di are measured on the lattice [27] and Di is found to be much 
smaller than D. Therefore one can neglect the nonperturbative part of V^{r), 
while that of V4 turns out to be also small numerically, Mss{nonpert.) < 
30MeV, and we shall also neglect it. 

The only sizable spin-dependent nonperturbative contribution comes from 
the term ^ (Thomas precession) and can be written at large distances as 

AH(Thomas) = -^I^ (25) 
r 2/i^ 

Now we come to the point of perturbative contributions to spin splittings. 
The simplest way to calculate those to the order 0{as) (and this procedure 
holds true for quarkonia) is to represent perturbative gluon exchanges by the 
same Eichten-Feinberg formulas (24) where one should keep in Vi,{r) only 
perturbative contributions to correlators D and Di in (A3.8)-(A3.11); then 
to the order 0{as) one obtains 

= 0, = ^ 2 ' 26 

r dr dr r^ 

^Apert) _ SC2{adj)as , . 

- ' l^'J 

= 87rC2iadj)as6^'\r) (28) 

However this procedure should be corrected for glueballs since i) valence 
and exchanged gluons are identical and ii) there is a 4-gluon vertex in ad- 
dition. The corresponding calculations have been done in [28], which show 
that corrections amount to the multiplication in (26) with the factor 3/4 and 
in (28) with the factor 5/8. 

With the account of these corrections the corresponding matrix elements 
in look like 

^SL -^^^(71) (29) 



j^ipert) ^ 57r(72(arfj) ^^^^(3)^^^^ 



3/i, 



^(pe.i) ^ ^M^(^(3sWnS(2)n - S^^^S^^))) (31) 
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iFvom (PPP one can see that Mss can be written as 

5a. 



Mss = ^\Rm' (32) 



4/i, 



To make simple estimates, we shall neglect first the interaction due to 
PGE between valent gluons. Indeed we show in the next section that this 
interaction cannot be written as Coulomb potential between adjoint charges, 
and comparison to perturbative BFKL Pomeron theory [12] shows that it is 
much weaker than Coulomb potential. Neglecting this interaction altogether, 
one gets the lower bound of spin-dependent effects, since all matrix elements, 
like {5'^^\r)), {-), {^) are enhanced by attractive Coulomb interaction. 

For purely linear potential one has simple relation, not depending on 
radial quantum number n 



Att Atx Air ^ ^ 

Using (33) and Mq = A^q, and taking Mq from Table 1, one obtains 

Mss = (34) 

and for rir = 0,1 and as = 0.3 the spin-spin splitting is 

Mssirir = 0) = O.SGeV, M,,(n^ = 1) = 0.20GeV (35) 

For M{0^~^) and M(2"^+) one has the values given in Table 5 for af = 
O.lSGeV"^ and for the sake of comparison with lattice calculations in Table 6 
for af = 0.228GeV^ and a, = 0.3. 

For L > one needs to compute spin corrections Msl and Mt- First of 
all one can simplify matter using the equation (it is derived in the same way, 
as (33) was derived in [29], for details see Appendix 4)) 



L{L + l){^,) = ^{V'{r)) (36) 



For V{r) = Uadj^ both M^'^'^^^ and m)?^''*'' are easily calculated and listed in 
Table 7. 

The nonperturbative part of spin splittings is due to the Thomas term, 
(^ + calculated numerically using the exponential form of D, Di 

found on the lattice [27], for details see [25]. 
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The resulting figures for AMThomas are given in Table 7. Combining 
all corrections and values of Mq from Table 2,3 one obtains glueball masses 
shown in Table 5 for a/ = O.lSGeV^ and compared with lattice data in Table 
6 for af = 0.228G'e1/2 

One can see in Table 6 that calculated spin sphttings of lowest levels 
are in good agreement with lattice data. This is another phenomenological 
manifestation of the PGE suppression in the glueball system; indeed had we 
taken PGE in the adjoint Coulomb form with ag = 0.3, we would obtain 3 
times larger spin splittings [18]. 

The general feature of spin-dependent contribution AHg is that it dies 
out very fast with the growing orbital or radial number, which can be seen 
in the appearance of the Hq factor in the denominator of (29-31). 

Indeed, from (8) one can derive that Mq ^ 4/io and therefore AHs ~ 
M2^rL) (^(r))' '^here O stands for terms like const. ^ or const'. ^ (from per- 
turbation theory). Hence spin splittings of the radial recurrence of states 
0+"*", 2++ or 0~+, 2~+ should be smaller than the corresponding ground states. 
This feature is also well supported by the lattice data in Table 5. 



4 Perturbative gluon ladders and glueballs 

In many analytic calculations of glueball masses it is postulated that there is 
a Coulomb-type interaction between valence gluons, which differs from the 
qq case by the Casimir factor, C2{adj) — 3 instead of C2{fund) = |. Before 
going into the details of the question how the perturbative gluon exchanges 
give rise to the Coulomb kernel, we here first assume that this is indeed the 
case, and correspondingly calculate the eigenvalues of the hamiltonian 

H^H^_C2{a^ (37) 

where Hq is given in (8). The resulting masses are hsted in Table 8 (the first 
three lines) for a, = 0, 0.2, 0.3, 0.39. 

One can see a drastic decrease of the mass due to the Coulomb attraction, 
especially for L=0. For a conservative value ag = 0.3 this mass drops down 
by 0.5 GeV. 

This is much larger than in the qq case [13,19], evidently due to the large 
Casimir factor. 
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Another characteristics of the Coulomb shift, which is useful for the com- 
parison with perturbative Pomeron approach [12], is the Regge slope 0^(0) 
and Regge intercept adO) of the glueball trajectory drawn as a straight line 
through the glueballs with L = 0(2++) and L = 0(4++) 0. These values 
are given in the last two lines of Table 8, and show a drastic increase of the 
intercept due to Coulomb interaction by AadO) ~ 0.64 for = 0.3. 

This will be compared later in this chapter with a similar large shift of 
the perturbative Pomeron trajectory Aap(O) in the lowest 0{as) approxi- 
mation [12] and with much smaller value of Aap(O) in the next (one-loop) 
approximation [30]. 

This comparison casts one more doubt on the validity of the assumption 
about the presence of the adjoint Coulomb interaction in the form (37). 

A similar conclusion can be deduced from spin-averaged eigenvalues. In- 
deed for the Hamiltonian 

, C2{adj)as 

n = Hq 

r 

the eigenvalues are given in Table 8. 

One can see that for L = both = 0.3 and 0.39 strongly contradict 
data, which shows that perturbative gluon ladder strongly differs from the 
adjoint Coulomb interaction, moreover the overall agreement of our results 
for Mo (where no Coulomb interaction is present) with spin averaged lattice 
masses tells that PGE is strongly reduced on the lattice. 

To study this point in detail one should consider the set of perturbative 
gluon exchanges and compare them to the BFKL diagrams describing the 
perturbative Pomeron [12]. 

First of all one should look into the mechanism which produces color 
Coulomb interaction, and it is instructive to compare quark-antiquark and 
gluon - gluon system from this point of view. For both systems there are 
diagrams of gluon exchanges in the order 0{g^), and in addition for gluon - 
gluon system there is in the same order the diagram of contact interaction, 
which affects the hyperfine splitting [28]. 

The main point is whether and how these diagrams are summed up to 
produce the color Coulomb kernel in the exponent, entering the Green's func- 
tion of the system. For the qq system (neglecting spin degrees of freedom for 

^This discussion is rather qualitative. Indeed Coulomb interaction modifies linearity of 
nonperturbative glueball trajectories. 
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(38) 



simplicity of comparison) one has the exact Feynman-Schwinger representa- 
tion 

Gqq^ j dsdsDzDze-^-^{W{C^^)) (39) 

where the Wilson loop is along the paths z, z integrated in (39). One can use 
the cluster expansion for purely perturbative gluons in W{C) as was done 
e.g. in [13,21] 

{W{C)) = exp[-C,^ / / + Dig')] (40) 

c c ^ ' 

For straight-hne trajectories z{t), z{t') (e.g. for static quarks) the integral 
in the exponent of (41) readily yields the color Coulomb potential, (W) ~ 

exp{^t). 

For light quarks one can consider the integral in the exponent of (40) as 
the full-fledged relativistic Coulomb kernel. It is legitimate to keep this kernel 
which is 0{g'^) in the exponent, neglecting Olg"^) additional terms provided 
the Coulomb kernel yields some amplification. 

This is indeed true in the nonrelativistic region (where Coulomb correc- 
tion are of the order -C 1) or at small distances (high energies) where 
this kernel yields double logarithmic terms [31]. Let us turn now to the gg 
system (the same is true a fortiori for the three-gluon system) . 

In (3) we have derived the gg Green's function for valent perturbative glu- 
ons in the nonperturbative background. The similarity of forms (3) and (39) 
is only superficial, since the Wilson loop in (39) contains both perturbative 
and nonperturbative contributions and one may argue, that perturbative ex- 
changes dominate at small distances, and hence exponentiate as in (40) and 
consequently give forth the color Coulomb kernel. 

In contrast to that , in (3) (Wp) contains only nonperturbative fields B^, 
yielding confining string between gluons, but no perturbative exchanges at 
all. In the framework of the background perturbation theory the perturbative 
vertices O(a^) and O(a^) enter the interaction Lagrangian, and there is a priori 
no guarantee that gluon exchanges produced by these vertices exponentiate 
to give a color Coulomb kernel. (Note that there is a difference between 
gluon exchanges, and spin - dependent vertices considered in the previous 
chapter, since the latter are taken as a perturbation in the lowest order, and 
there is no need for them to exponentiate to the Coulomb ladder). 

Having all this in mind, we turn our attention to the subset of graphs 
which is summed up in the BFKL approach [12], and which is distinguished 



14 



by the principle of leading diagrams in the high-energy scattering, or in 
another setting, by the summation of ladders for the leading Regge trajectory 
in the t-channel. Since these ladders are dominant perturbative series (see 
[12]) for the Pomeron trajectory, we can consider the same contribution in 
our circumstances - for calculation of glueball masses, extending in this way 
BFKL - type analysis from Pomeron-generating glueballs (4++, 2++eic) to 
all others, and having in mind, that this may give only an estimate of the 
order of magnitude. 

Thus our purpose is now to estimate the contribution of the BFKL dia- 
grams to the glueball masses (perturbative mass shift) and compare it to the 
usual color Coulomb contribution. 

In order to estimate effects of small distance contributions we shall use the 
analysis of these effects on gluonic Regge-trajectories not from the glueballs 
mass spectra at positive t, but for t = 0. Extensive calculations of the 
gluonic Pomeron trajectory intercept have been carried out in the leading 
log approximation (LLA) [12] and ag corrections were calculated recently 
[30]. It has been shown that the leading Regge singularity corresponds to a 
sum of ladder type diagrams, where exchanged gluons are reggeized. In the 
leading approximation an intercept of this singularity is equal to [12] 

47V 

ap{0)^l + as — -ln2 (41) 

71 

The shift from the noninteracting gluons point ap{0) = 1 is equal to A = 
ap{0) — 1 0.5 for ~ 0.2. This rather large shift is strongly reduced by 
as corrections [30] 

A = as—ln2{l - Ca,) (42) 

TT 

The coefficient C is rather large 6.5) and the as correction strongly re- 
duces A. Its value depends on the renormalization scheme and scale for a^. In 
the "physical" (ELM) scheme values of A are in the region 0.15-=-0.17 [30]. In 
this approximation the leading gluonic singularity is Regge-pole. and we can 
estimate mass-shift of the lowest glueball state using this result and assuming 
that the slope a'p = ^^-^ — ^ 0.46*6^^^(11) will not be strongly modified by 
perturbative effects. Thus one can expect that characteristic shift due to per- 
turbative effects in M\L = 0, rir =_0), 6M^ ^ ^/a'adj ~ (p-38 ^ 0^48)^6^^ 
This corresponds to the shift in M{L = O^rir = 0),6M f» SM'^/2M 
O.lGeV. This shift should be compared with much larger mass shift from 
pure Coulomb interaction given in Table 8. Thus the 0{as) correction to 
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the BFKL ladder gives a strong suppression of PGE series and may be a 
possible explanation why Coulomb -like attraction is not seen neither in 
spin-averaged masses M{L,nr), nor in spin splittings. It should be noted 
that this is only a rough estimate of the perturbative effects because higher 
orders of the perturbation theory can modify this result. 



5 Three— gluon glueballs 

The three-gluon system can be considered in the same way, as it was done 
for the two-gluon glueballs. The 3g Green's function G*-^^^ is obtained as 
the background-averaged product of 3 one-gluon Green's function, in full 
analogy with (1). Assuming large Nc limit for simplicity and neglecting 
spin splittings and projection operators one arrives at the path integral (cf 
equation (3)) 

G^^^) = const TT / dsiDz^^e-""'-"^' (43) 

where a = a fund, since every gluon is connected by a fundamental string 
with each of his neighbors. 

Using as before the method of [13,14] and three-body treatment of [17] 
one obtains omitting spin-dependent terms the following Hamiltonian (we 
assume symmetric solution with equal //^(t) = //(t),^ = 1,2,3 (no orbital 
excitations was assumed as in (8)) 

a,3»,^P?±P| + | + , (44) 

^/^ ^ i<0=l 

Here r^j = |rj — r^ l, and is the space coordinate of the i-th gluon, while 
^, p^i and 77, are defined as 

ri-r2 /3ri+r2 Id Id 

To simplify treatment further, we shall consider ^ clS cl constant to be 
found from the extremum of eigenvalues, as in (8), which in that case provided 
some 5% increase in eigenvalues (see Table 3 of [19]), and what we expect 
also in this case. 
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To find eigenvalues of H^^^^ one can use the hyperspherical method in- 
troduced in [32] and appUed to the 3q system in [17]. Defining hyperradius 
p, = 77^ + ^ , one obtains one-dimensional equation for the eigenfunc- 
tion XniP) = x{p)j{K is the grand angular momentum K = 0, 1, 2..., and n 
-radial quantum number). 

-^X" + Ueff{p)x{p)-Mx{p) (46) 



where 



,2 , . 15^ , 32^2 



Solution of (46) is expressed through generalized Airy functions. 

A reliable (within few percent of accuracy) estimate of M is obtained 
when one replaces Ucff{p) by the oscillator well, with the center at the min- 
imum of Ueff{p), p = Po, and frequency luq expressed through U"j-jr{po). For 
K — one has 

In this way one obtains 

and minimization of M[p) in p yields 

PO = {W\^Y"; M{p,) = 6/.0 + f = 6/.0 + (^)3/^(^) V2 (50) 

For a = O.lSGeV^ one obtains loq = l.lSGeV and po = O.UGeV, hence 
the minimal eigenvalue is 

Mo = 3.23GeV (51) 

This spin-averaged value is listed in Table 4. Radial excitations are given by 
an approximate equation 

= Qpo + Y + ncoo, n = 0,1,2... (52) 



Orbital excitations yield increase in mass of the order 

AM-^^;±i^, (53) 
2ppl ^ > 

17 



which for lowest excitation gives 



AM^=^ ~ IGeV (54) 

almost the same amount, as for the radial excitation. 

Coulomb shift (if Coulomb interaction existed between gluons) , would be 
enormous: AM{Coulomb) = —1.3GeV. However here one can use the same 
arguments as for two-gluon glueballs and drop the color Coulomb interaction 
between gluons altogether. 

Finally we turn to the question of quantum numbers and spin splittings 
of the 3g states. 

According to (27), perturbative hyperfine interaction is given by matrix 
elements 

AM,, = E(S«S(^)^^^^^a.5(3)(r.,)) (55) 
i>j -^/^o 

Note that in the large Nc hmit gluon lines are replaced by double fun- 
damental lines and planar gluon exchanges occur with fundamental charge, 
hence the fundamental Casimir operator in (55). 

For the K = state the wave function depends only on hyperradius p 
and one has 

(^^'Hr.)) = (^) (56) 

Now for K = state all internal angular momenta are zero and one can 
express (SjSj) through total angular momentum J 



(S.S,) = iil+ltll (57) 



As a result one obtains AM55 for — 0.3 and a — O.lSGeV'^ 



AMss = — '^,3 ^ ' (S,S,) ^ 0.644Ato^ ^ GeV (58) 



Hence the candidate for the odderon state J^^ = 3 is shifted by 
0.286*6^^ upwards and 1 by 0.189G'ey downwards with respect to (51). 
Resulting values of glueball masses are listed in Tables 5 and 6. 
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6 Glueball Regge trajectories and Pomeron. 



The leading Regge trajectory (with the largest intercept Q;p(0)) is usually 
called the Pomeron trajectory. It plays a special role in the reggeon approach 
to high-energy hadronic interactions. The parameters of the Pomeron tra- 
jectory and especially its intercept play a fundamental role for asymptotic 
behaviour of diffractive processes. We already touched the problem of the 
Pomeron intercept in Section 4, where the perturbative, small distance con- 
tribution has been discussed. Now we will consider this problem in more 
details, taking into account both nonperturbative and perturbative contri- 
butions to the Pomeron dynamics. 

The large distance, nonperturbative contribution gives according to Eq.(12) 
for the leading glueball trajectory (n^ = 0) 

ap(t) ^ -ci + a'pt (59) 



Taking into account spins of "constituent" gluons but neglecting small 
nonperturbative spin effects we obtain for the intercept of the leading tra- 
jectory 

ap{0) = -ci + 2 (60) 

which leads to ap{Q) 0.5, and this value is substantially below the value 
found from analysis of high-energy interactions ap{0) = 1.1 ^ 1.2 [20]. 

The perturbative (BFKL) contribution leads to a shift (increase) of the 
Pomeron intercept by 0.2, as it was explained in Section 4. The resulting 
ap\o) 0.7 is still far from experiment. 

There are other nonperturbative sources, which can lead to an increase 
of the Pomeron intercept. One of the most important is in our opinion the 
quark-gluon mixing or account of quark- loops in the gluon "medium" . In the 
1/N -expansion the effect is proportional to Nf/N^., where Nf is the number 
of light flavours and this mixing is known to be important (at least in the 
small-t region). 

In the leading approximation of the 1/iVc-expansion there are 3 Regge- 
trajectories, - gg-planar trajectories {f^^^ made of uu and dd quarks and 
/'^°) made of ss-quarks) and pure gluonic trajectory - G. The transitions 
between quarks and gluons ~ will lead to mixing of all these trajectories. 
Lacking calculation of these effects in QCD we will consider them in a semi- 
phenomenological manner. From the mixing of two trajectories 1 and 2 with 
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the transition constant gu it is easy to obtain the following values of new 
trajectories (see for example [33]) 



(i) ai{t) + a2{t) ± ^{ai{t) - a2{t)y + AgUt) 
«± = ^ (61) 

Note that for a realistic case of G,f and /'-trajectories (Fig.l) all 3- 
trajectories before mixing are close to each other in the small t region. Tra- 
jectory of gluonium crosses planar / and /'-trajectories in the positive t 
region {t < IGeV^). In this region mixing between trajectories is essential 
even for small coupling matrix gik{t)- 

The dual unitarization scheme [34-36] leads to the conclusion that the 
quantity gl2 fastly decreases as t increases in the positive t region. This 
means that at large positive t, a±{t) coincide with trajectories ai and a2, as 
it happens in (61) for gl2 <^ |q;i — q;2| 

2 2 

a+^ai+ ; a- a2 H (62) 

ai — a2 a.2 — cti 

This phenomenon is called the asymptotic planarity [36]. Let us note 
that mixing effects will be small in the large t region even if the couplings 
have weak i-dependence because the differences between planar and gluonic 
trajectories increase ~ i at large t. 

For weak mixing between trajectories {g^^ <^ \af^ — oc^^\) (62) can be 
generalized 

^ «f + E (0)^^' (0) (63) 

For g'j ~ 0.1 typical resulting trajectories are shown in Fig.l by solid 
lines. For details see appendix 5. The Pomeron trajectory is shifted to the 
values Q;p(0) > 1. For t > IGeV"^ Pomeron trajectory is very close to the 
planar / trajectory. 

Position of the second vacuum trajectory at t < is close to cij, while 
for t > IGeV'^ it is close to ctj/. The third vacuum trajectory is below afi 
at t < and at t > IGeV it is close to the ac- Due to asymptotic planarity 
effects of mixing are not very essential for properties of physical particles 
on these trajectories as all resonances are in the region t > l.bGeV"^. On 
the other hand they are important for understanding of the SU(3)-breaking 
effects for the Pomeron exchange amplitudes at t < 0. 
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At the end of this section let us consider the "odderon"- the leading 
Regge trajectory with negative signature and C-parity. Mass of the lowest 3g 
glueball with spin 3 corresponding to this trajectory has been estimated in 
the previous section and found to be large ~ 3.51GeV ( for cx/ = O-lSGe^^) 
in accord with lattice data. The slope a'^g for this trajectory should be equal 
to the one of (7(7-trajectory0 and thus the intercept of the nonperturbative 
glueball "odderon" is very low Q!3g(0) ~ —1.5. Mixing with gg-trajectories 
{uj, Lf) is much smaller than in the Pomeron case since there is no crossing 
of the odderon and (cj, Lp) trajectories in the small t-region and thus gluonic 
"odderon" is not important for high - energy phenomenology (at least in the 
small-t region). 

7 Discussion and conclusions. 

The main results of the paper can be separated in two groups. In the first 
part we calculate the 2g and 3(7 glueball spectrum analytically and compare 
resulting masses with lattice data, finding a very good agreement. In the 
second part the glueball Regge trajectories are obtained and their correspon- 
dence with the Pomeron and odderon is discussed. 

Concerning the glueball spectrum the spin-averaged results of section 3 
calculated for all states of 2g and ?>g glueballs yield a very good agreement 
between our results and spin-averaged lattice masses. We stress that our 
spectrum in contrast to most existing theoretical models contains no fitting 
parameters, and all masses are expressed in terms of string tension a, as it 
is done also on the lattice. 

This coincidence and evident smallness of PGE interaction, which would 
be very strong if it had been adjoint Coulomb interaction due to the Casimir 
factor 3, has called us for more detailed investigation, whether Coulomb is 
indeed appropriate in the systems of valence gluons. The analysis done in 
Section 4 has brought us to the conclusion that the situation in the system 
of valence gluons is completely different from that of valence quarks, and the 
perturbative gluon exchanges do not exponentiate into Coulomb kernel for 
2g and ?>g systems in contrast to the gg, 3g systems. 

This observation explains qualitatively the absence of strong Coulomb 

•^The situation is analogous here to the case of qq (meson) and qqq (baryon) Regge 
trajectories, where baryon trajectory displays the quark-diquark structure and hence the 
meson Regge slope [17]. 
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downward shifts of glueball masses and moderate spin splittings in lattice 
calculations. To make a quantitative estimate, we have considered the BFKL 
perturbative series for Pomeron [12] including one loop correction [30]. This 
series is not a Coulomb ladder and with an account of NLO corrections it 
leads to a mass shift of approximately 3-4 times smaller than for Coulomb 
interaction. 

In contrast to Coulomb interaction the spin splittings of glueball masses 
are obtained from the first perturbative correction calculated with nonpertur- 
bative wave functions. There is a good agreement on spin splittings (within 
few tens of MeV) between our calculations and lattice data, as it is shown in 
Table.6. 

The agreement signifies that the main ingredient of the glueball dynamics 
is the adjoint string (or two fundamental strings) occurring between gluons in 
the two-gluon glueballs and the triangle construction of fundamental strings 
in the 35^ glueballs. The string dynamics implies that glueball masses lie 
on the corresponding straight-line Regge trajectories, which have the Regge 
slope equal to | of that for meson trajectories. In other respects glueball 
trajectories are similar to the qq trajectories for zero-mass quarks [13-16, 
19,24]: they are straight with good accuracy, have daughter recurrences due 
to radial excitations, which are also approximately straight lines. 

In the last part of the paper we used our knowledge of glueball Regge tra- 
jectories for investigation of the Pomeron singularity. The Pomeron, which 
yields asymptotically dominant contribution at large energies, is certainly a 
complicated object, which has some features connecting it with the dominant 
glueball trajectory. First of all, Pomeron exchange has a cylinder topology 
(which is supported by the multiplicity analysis [37]), similar to that of glue- 
ball amplitude which becomes evident when one replaces the adjoint string 
by the double fundamental string. 

The idea of Pomeron as a two-gluon exchange amplitude has a long his- 
tory [11] and was exploited both in perturbative [12], nonperturbative [18] 
and hybrid [38] approaches. The purely perturbative approach has some 
difficulties of internal consistency both because of slow convergence of per- 
turbative series [30] and because of sensitivity to large-distance contributions 
[39]. The latter signifies that the nonperturbative effects may play very im- 
portant role in the Pomeron dynamics, and our paper is a demonstration 
of this. The character of nonperturbative trajectories is linear due to the 
string dynamics and absence of mass dimension parameter other than string 
tension. Perturbative singularities in the j-plane are not always poles, and 
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are certainly not linear trajectories. 

Our discussion in previous sections has come to the conclusion that per- 
turbative effects only slightly shift nonperturbative trajectories (increasing 
Regge intercept by~ 0.2). 

With all that we have come to the Regge intercept of leading Regge trajec- 
tory equal to 0.7, strongly different from the experimental Pomeron intercept 
of 1.07 ^ 1.2. And here comes an interesting observation, made earlier in a 
bit different context [33], that due to different slopes of meson and glueball 
trajectories they must intersect each other in the region of small t < IGeV^. 
We have taken this fact into account in the three-pole model, where coupling 
constants between channels /, /' and G are introduced phenomenologically. 

The results, shown in Fig.l demonstrate a dramatic change in the course 
of trajectories: the largest intercept increases by 0.5 reaching the physically 
reasonable values of 1.2 

Let us note that both nonperturbative (string dynamics, quark loops) 
and perturbative effects are important to obtain ap{0) > 1. It is impossi- 
ble to separate out "soft" and "hard" Pomerons, as sometimes is done in 
phenomenological studies of high-energy interactions of hadrons and small-x 
physics of deep inelastic scattering. 

In a supercritical Pomeron theory with A = q;p(0) — 1 > the cor- 
responding multipomeron exchanges are important at very high energies. 
They allow to obtain scattering amplitudes, which satisfy s-channel unitar- 
ity and the Froissart bound for total interaction cross sections as s — oo. 
From the point of view of 1/Nc - expansion multipomeron exchanges are 
~ {l/N^Y^ , where n is the number of exchanged Pomerons, but they have 
faster increase with energy (~ (~)"^) than the pole term and should be re- 
summed. This can be done using Gribov's reggeon-diagram technique [40]. 
In practical applications of Reggeon theory to a description of high-energy 
hadronic interactions multipomeron exchanges are essential for simultaneous 
description of total interaction cross sections and multiparticle production 
(for review see [41]). 

Looking back to a structure of vacuum trajectories we found that each of 
3 new trajectories ai{t) is now a mixture of G, /, /' and only asymptotically 
at large t they tend to the original trajectories. As it is seen in Fig.l the 

^Note that we discuss here the "bare Pomeron" intercept, which is larger than an 
"effective" Pomeron intercept usually extracted from HE scattering. As was discussed in 
[20] the bare Pomeron characteristics are measured in small-x DIS experiments, yielding 
the intercept around 1.2. 
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leading trajectory (with largest intercept), which should be associated with 
Pomeron, asymptotically tends to /, the second trajectory as positive t is 
close to /', while the third asymptotically (at large t) coincides with G while 
at t = it is below the first two trajectories. Thus a rearrangement takes 
place: G trajectory is shifted downwards, while the / trajectory is lifted up 
and becomes the Pomeron. 

One of immediate consequences of this rearrangement is the special pat- 
tern of Pomeron couplings, which can be measured experimentally. While G 
trajectory was flavour blind, now due to mixing one can calculate the cou- 
plings of the Pomeron to light quarks (via /) , to strange quarks (via /') and 
symmetrically to all flavours (via G). 

Mixing between gluons and qq pairs has another important aspect, - it 
leads not only to shifts of Reai{t) but also to appearance of Jmacit) and 
as a consequence to nonzero widths of resonances on glueball trajectories. 
They should be of the same size as mass shift due to the mixing and thus 
are expected to be not too large, Pg ~ lOOMeV. 

Present study can be improved in several points. Firstly, perturba- 
tive contributions to the glueball trajectory including spin-dependent terms 
should be studied more systematically. 

Secondly, analytic calculations of gik{t) are necessary to make our theory 
complete. Finally a detailed analysis of experimental consequences of our 
results is needed, which is planned for separate publication. 

The work of A.B.Kaidalov was partly supported by the grants RFBR 
98-02-17463 and NATO grant OUTR.LG 971390. 
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Appendix 1 

Creation operators of glueball states. 



We consider in this Appendix and Tables 9,10 the operators vl/^'"-^ and 
^(owi) (A2.16)-(A2.17) which specify glueball states and their quan- 

tum numbers, J^^ . One may consider local \t'(x, x) or nonlocal operators 
^(x, y) for two-gluon glueballs and corresponding operators for many-gluon 
glueballs '^{x^^\ ...x^")). For simplicity listed below are only local versions, 
since nonlocal ones can be constructed with the help of parallel transporters 
^{x,y), as it is done in (A2.16), (A2.17). 

First one can construct ^'^ in a general form, not assuming separation 
of into background and valence parts, similarly to what is done on the 
lattice. Then one has at his disposal vectors Ea, Da pseudovector Ba and 
color tensors 5ab, fate, dabc- One should also take into account that under 
charge conjugation C the following transformations hold 

CF,,C-' = -Fl, CD,C-' = -Dl (Al.l) 

Hence one obtains the following list of states for the two-gluon glueballs 
(containing two field operators) and due to Bose statistics symmetric with 
respect to exchange of all coordinates of two gluons. We also hst in the first 
column of Tables 9,10 the dimension of the corresponding operator. 

In the background perturbation theory (BPT) and vf^^""*) can be 

constructed from the spacial components of the gluonic field ai,i = 1,2,3 
since the fourth component 04 can be expressed via the background gauge 
condition 2?^a^ = 0. Note that Oj transforms homogeneously (see equation 
(A2.4) of Appendix 2) and therefore one obtains gauge-invariant combina- 
tions for vI>(o"*) replacing Ei in the third column of the Table 9,10 by 
whereas J^^ does not change. In the same way is replaced by (D x a)^, 
and one obtains the fourth column of the Table 9,10. Dimension of BPT 
operators is given in column 5 and the orbital angular momentum in the last 
column. 

For three-gluon glueballs the corresponding entries are given in Table 10. 
One should notice that C parity of all listed states is here negative. Again, 
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dimension of BPT operators is listed in the last column. 

As can be seen from the results of our calculations in Tables 4-6, the 
glueball spectrum is in good agreement with the hierarchy associated with 
increasing angular momentum L or BPT dimension (they differ by two units 
for two-gluon glueballs) The same ordering persists in lattice data. Three- 
gluon glueball masses are typically shifted by 1.5 -^ 2GeV (an exception of 
1"'"~,3^~ states in lattice data waits for explanation). Note the absence of 
J = states in the lattice spectrum. From Table 9 one can see the only 
candidate, l""*", but the corresponding local operator is proportional to the 
energy-momentum tensor and by the arguments of [42] the residue of the 
state should vanish. 

In both Tables 9,10 we have used following notations 

symrrHik)Tik = Tik + Tki - '^SikTu 

Symmetrization of higher operators TikiiTikim is done in the usual way to 
construct irreducible 0(3) tensors. 

Appendix 2 

Glueball Green's function and Hamiltoman in the background 
formalism. 



In what follows the Euclidean space-time is used. 

The total gluonic field ^4^ is spht into NP background and valence 
(perturbative) gluon field a^, 

A^ = B^ + a^. (A2.1) 
The QCD partition function Z{J), 

Z{J) ^^J e-^^(^)+/'^^(^)^^(^)'^'^D^L>V'£>V' (A2.2) 
where Se is can be rewritten using 'tHooft identity as 

Z{J) = ^ / DBr]{B)e! ^'''^^^ ( DaDx/jD^j^e-^^^^+'^'^+f '^"^''^ (A2.3) 
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Here r]{B) is (an arbitrary) weight of integration over background fields 
B^, the exact form of which is not of interest to us, since the overall effect 
of background fields will enter in our results via string tension a and (in 
some corrections) as a NP field correlator {F{x)F{y)). Both quantities are 
considered as an input. 

In what follows we shall expand (A2.4) in powers of ga^ as it is usually 
done in background perturbation theory [21,43]. In the lowest order of l/N^ 
expansion quarks are decoupled from gluons and we shall neglect coupling 
to quarks till the last two sections of the paper. 

It is convenient to prescribe the following gauge transformations, 

^ U+a^U, (A2.4) 

^ U{B^ + '-d^)U, (A2.5) 

and to impose on background gauge condition 

D,a, = d,a1 + gr^'By^ = (A2.6) 

In this case the ghost fields have to be introduced and one can write the 
resulting partition function as 

Z{J) = ^ / DBr]{B)exp{ J J^B^d^x)Z{J, B) (A2.7) 



where 



Dadet{jj)exp J d^x[L(a) - -{G^f + J^a^] (A2.8) 



where 



L{a) = Lo + Li{a) + L2{a) + Lint{a) 
L2{a) = +^au{DlSi^i, - D^Du + igF^y)a^ = 

^atiD'^^DlX, - D^;D':' - gr'F^Mi (A2.9) 
D^^dx-5ca + gr''''B'^ = bx 
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Lint = -\{D,{B)a, - D,{B)a,rgr'%al - -/ r'^air^ alal (A2.10) 
The background gauge condition is written as 

= d,al + gr^'Blal = {D^a.f (A2.11) 
and the ghost vertex [21,43] is obtained from ^^(Di^(B)Dn(B + a))ab to be 
Lghost = -Qt{D^{B)D^{B + a))abQb (A2.12) 

The hnear part of the Lagrangian Li disappears if B^ satisfies classical 
equations of motion. 

We now can identify the propagator of from the quadratic terms in 
Lagrangian ^2(0) - ^(G")^ 

Gt = [DlS^.. - b,D. + igF^. + -^D.D^a^ (A2.13) 

It will be convenient sometimes to choose = 1 and end up with the 
well known form of propagator in - what one would call - the background 
Feynman gauge 

Gfu = {Dl ■ 5,. + 2igF^,)-' (A2.14) 

We are interested in the glueball Green's function and therefore must 
define first the initial and final state vectors of glueballs, consisting of 
valence gluons initially and oi nj gluons in the final state. The following 
general nonlocal state vectors can be used for k gluons 

...x^^-^^) = r[/o(a(x(°))0(x(°),x«)/i(a(x«), .../fe_i(a(x('=-^))(/.(a;(^-i), 

(A2.15) 

X 

Here (f){x,y) — P exp (ig J Bn(z)dzij) is parallel transporter, all are 

y 

in fundamental representation, and /(a) is a polynomial in a^, which may 
contain derivatives in the form of = 5^ — igB^. 

According to (A2.4),(A2.5) are color singlets. One can also have local 
form of ^'fe taking all a;'-*-' at one point. Exact form of ^'^ is given in Appendix 
1. 

As will be seen below the state (A2.15) will evolve as a closed fundamental 
string with k gluons "sitting" on the string, when all /j are linear (and more 
gluons, when some /j have larger power). This form of initial and final 
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states is convenient for multigluon glueballs and will be used for three- gluon 
glueballs in Section 5. 

Another form of \E'fc (equivalent to the previous in the limit Nc — > oo) 
obtains when one takes adjoint string. E.g. for two-gluon glueballs one has 

^2(a;W,a;(2)) =i[a^(a;W0(a;W,a;(2)^)^^^^(2)| ^^2.16) 

Here the hat signs imply adjoint representation. For two-gluon glue- 
balls we are using (A2.17) for initial and final states, and the corresponding 
Green's function describes the evolution of the open adjoint string with ad- 
joint charges (gluons) at the ends. 

It can be calculated using (A2.4) in the form given by equation (1) of 
the main text, where we neglect terms Li{a) and Lint{a) (the first gives an 
insignificant correction discussed in [21], while Lint contains higher powers 
of (yfa^j, and will be used for calculation of perturbative corrections to the 
Green's function). 

The next step is the Feynman-Schwingcr representation (FSR) [13,14] 
for the gluon Green's function (A2.15), which allows to exponentiate and 
F/^i^ as 

X s 

y (ig J Bi,dz^+2ig J F{z{T))dT) 

G^^{x,y) = const dsDze-^PePpe ^ ° (A2.17) 

oo 

Insertion of (A2.17) into (1) and using the fact that ordering inversion for 
one of gluons yields (— -Bj, —F'^) instead of i?^, F brings forth the equations 
(3) and (4) of the main text. 

Equations (3), (4) is another form of (1),(2) and contain no approxima- 
tions (except for omitting of the terms Li and L^^t discussed above and used 
in writing (1) and (2)). 

Another important step done first in [13,16], and developed in [14,15], 
is the introduction of the auxiliary function /i(r), which may be called the 
einbein, and which plays a crucial role of effective gluon mass in the whole 
formalism. This is done rigorously and without introduction of arbitrary 
fitting parameters in contrast to usual potential models. 

Defining 

2/x(t) = = ^4 (A2.18) 
or 
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where r or s is the Schwingcr proper time, and t is the Euchdean time at 
any point of trajectory Zh{t), one can rewrite FSR (3) identically as 

G^,,^^,^{x,y\x',y') = const J Dii{t)Dii{t')DziD4e-''-''' {Wp) (A2.19) 

Where Dzi (or Dz'^) denote 3d path integrals over trajectories Zi{t),i — 1,2,3 
and z'i^(f/), k = 1, 2, 3, and D^(t) is the Id path integral over functions //(t). 
Kinetic terms K, K' can be expressed through e.g. 

T 

J^m{t))' + mT^x,-y, (A2.20) 



where dot means time derivative. The form (A2.20) reminds nonrclativistic 
kinetic energy however is exact relativistic form. In case of massive rela- 
tivistic particle with mass m in the corresponding term in action, K^, looks 
hke ^ 

Km = /(^ + ^amr + 1])^^ (A2.21) 

Introducing momentum Pj = §^^, one would obtain after extremization 
in the usual answer for the Hamiltonian 

Ho = -^1+^ (A2.22) 

In case of zero mass, m — 0, one would obtain for the free gluon without 
spin from (A2.21) the free hamiltonian Vp^- 

The NP interaction in the two-gluon system is given by (4), where the 
term J B^dz^j^ generates the adjoint string, equation (5), and after introduc- 
tion of another einbein function as it is done in [14] one obtains the 
Hamiltonian (6). 

The latter describes the straight-line adjoint string connecting two gluons, 
which can rotate and change its length. 

Contribution of NP spin terms is considered in Appendix 3. 



Appendix 3 
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Nonperturbative spin— splitting terms 



Introducing the spin matrix of the gluon as in (22) and using (A2.18), 
one can rewrite the F terms in the exponent of (4) as 

where (SF) = SiB^ + SiEi = S^^Fi^^. Note that the contribution of the 
second (primed) gluon to (A3.1) has another sign as compared to the first 
gluon. This is the consequence of the fact, that color and time ordering of 
operators B, F and B', F' are opposite in the closed loop Wp- One must use 
therefore in Wp the transposed operators for the second (or the first) gluon 
and write B'^ ^ -B' . 

To calculate the average of the exponential (4) , one can use the following 
trick: using the nonabelian Stokes theorem, and cluster expansion in the 
Gaussian approximation one rewrites Wilson loop integral as 

{Wadj) = {Ppexpig j F^^{u)da^^{u)) = exp[-^ J J {F^^Fxa)da^^daxa\ 

s s s 

(A3.2) 

Then (4) can be rewritten as 

T T 

iW,) = tre.pl-2. / ^S..^ + 2, J 

(A3.3) 

Evaluating derivatives, one arrives at the expression, based on the Gaus- 
sian approximation 

r r dt dt' 

{WF)=trexp[-^ J y[(ia^,(«)-2^5^,^][daA.(«')+2^^L^^](^M-(«)^Aa( 

(A3.4) 

One can add in the exponent of (A3. 2) all higher correlators in spin-independent 
terms and thus restore the area law (5) with aadj exact (i.e. beyond Gaussian 
approximation). For spin-dependent terms higher correlators bring about 
higher powers of S, S'. 

Since spin-dependent terms are relatively small corrections, it is legiti- 
mate to keep for them lowest, i.e. Gaussian approximation and write 

(Wf) ~ Ztrexp{-aadjSmin)exp{Ni + N2 + N^) (A3.5) 
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where notations used are 

N, = ig^ j J daxM ■ S^,^^{F^MF>.aHt))) (A3.6) 

and N2 is obtained from A^i by an exchange t ^ t' . 

The transformations in (A3. 6), (A3. 7) into the spin-orbit, spin-spin and 
tensor terms in (24) are the same as in the corresponding heavy quarkonium 
expressions given in [25], which are similar to (A3. 6), (A3. 7) modulo numerical 
coefficients and different gluon spin factors. 

The field correlators {FF) enter the final expression via potentials Vi{r),i — 
1,2,3,4 which are the same as for heavy quarkonia and given in [25], with 
the replacement C2(/) — > C2{adj). If one introduces two scalar functions D 
and Di as in [23] , one can write, 

IdVi 7 , fdX,^ A 



- / dvJ-{l--)D{X,v), (A3.8) 



R dR 

-00 



R dR 



-00 

00 

V3 = -/,„«^5^M|£) (A3.10) 

—00 

7 dD 
Vi= J dv[3D{R,v) + 3D^{R,v) + R^-^] (A3.11) 

—00 

Note that D = ^^D^, and the same relation for Di, where , D( 
refer to the fundamental representation. The normalization of D can be 
obtained from the relation 

^^aIj-\jd'xD{x) (A3.12) 

where the subscript (2) in (Jadj denotes the lowest (quadratic) correlator 
contribution to the string tension. As one can argue, the accuracy of this 
quadratic approximation is around 10% [23]. 
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Taking asymptotically large R in (A3. 8), and using (A3. 12) one obtains 

the asymptotics (25), given in the main text. 

To evaluate NP spin-orbit splitting, we must estimate the matrix element 
{V( + V2), i.e. some integrals with D, Di. The latter have been measured on 
the lattice [27] and found to be of exponential form 

D{x) = D{Q)exp{-x6), Di{x) < D{x) (A3.13) 

with S IGeV. For an estimate of {V{ + V2) we neglect Di, and calculate 

{VI + V^) = ^(2[- Ji(x) + ^""^'^^^"^^ 1 + Mx)) (A3.14) 

TT X 

where Jn{x), Kn{x) are Bessel and McDonald functions respectively and x — 
6r. 

^Prom (A3. 14) one can see that asymptotic behaviour (25) is obtained 
only for large r > 76"^ ~ 1.5/m. Therefore the average of AH{Thomas) 
with the square of the glueball eigenfunction is considerably reduced as com- 
pared to the average of the asymptotics (25), and the resulting NP spin- 
orbit term given in Table 7 is smaller than the corresponding perturbative 
term, and the ordering of the levels is due to the perturbative part of spin- 
dependent forces. 



Appendix 4 

Derivation of the relation for the matrix element (r~^). 

Writing solution of the Hamiltonian Hq, Eq.(8) in the form 

Mr) = ^Yi^, y„(r) ~ r'+\ r ^ (A4.1) 
r 

one can rewrite equation as 

y'; = [2fi{V{r) - E^) + ^^^^^]yn{r) (A4.2) 
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We shall use the procedure suggested in the second of ref. [29]. Multi- 
plying both sides of (A4.2) with and integrating over d^r, one obtains 

/ '^'3 = = JdrmV(r) - E„) + ^:<^l^f = 



= -I JyUr)drm'{r) - 2^^^^] (A4.2) 



Taking into account (A4.1) one has from two results: 
For L — one obtains the well-known relation [29] 

\MO)\'-^{V'{r)) (A4.3) 

For L > one has instead, 

L(L + l)(l)=/.(y'(r)) (A4.4) 

In our case (Eq.(8)) /jLq = 2/i, V'{r) = aadj- Note that both l.h.s in (A4.4) 
and (A4.5) do not depend on the radial quantum number n. 



Appendix 5 

Mixing of glueball and /, /' trajectories. 



One can start with the scattering amplitude of hadron a on hadron b in 
the Regge-pole approximation 

i,k 

where {i, k) = 1,2,3 refer to the bare Regge trajectories 

J = ai(i),i = l,2,3 (A5.2) 
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and the matrix Tj^ has the form 

T,, = {j - ait) - git))7k (A5.3) 

with the notations 

{a{t))ik = oii{t) ■ dik, {g{t))ik = gik{t),gii = (A5.4) 

The nondiagonal matrix g{t) describes mixing of Regge trajectories. In 
what follows we consider three bare trajectories: di{t) is the glueball trajec- 
tory calculated in the paper in Section 6. We approximate it in the region 
< t < 6GeV'^ by hnear form 

ai{t) = ai(0) + a[{Q) ■t^Q.7 + 0.246i (A5.5) 

where the value for a'i{0) is chosen in such a way as to reproduce the first 
glueball 2++ state at M = 2.3GeV (Table 5). 

The bare / and /' trajectories are denoted as 0:2 (i) and as{t) respectively 
and taken in the form 

a2{t) = 0.55 + 0.89t; a-i{t) = 0.25 + OMt (A5.6) 

The mixing matrix gik{t) is not known theoretically; as was discussed 
in section 6, condition of planarity [34-36] requires gik{t) to fall off at large 
positive t, therefore we assume for it the form 

(0) 

and for explicit calculations in the region t > Owe use k = 1 and = ^GeV^. 
To find shifted Regge poles in T, one can rewrite (A5.3) as 

Tik = TTTT (A5.8) 

det(j - a(t) - g(t)) 

where tik are minors of T. The roots of determinant in (A5.8) are given by 
cubic equation 

f-f^oii+j{^aiak-gik)-0£ia2a3+ J2 Qik^^i - '^912913923 ^ (A5.9) 

i^k ij^kj^l 
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We denote three roots of (A5.9) by 



j = ai{t),i = 1,2,3 



(A5.10) 



Let us start with t = 0. Assuming for aj(0) the values in (A5.5), (A5.6) 

(0) 



and the following values for g^' ^ 



{g?jr = 0.16; 



0.08; 



^13 



0.01 



we obtain the intercepts of mixed trajectories 

ai(0) = 1.2; asiO) = 0.075; a2(0) = 0.225 



(A5.11) 



(A5.12) 



Thus we obtain a realistic intercept of the Pomeron, corresponding to the 
bare Pomeron intercept observed in DIS' at small x [20]. Note however, that 
theoretical uncertainty in gl'^^ and in the Pomeron intercept are rather larger 
(~ 0.1). 

Resulting trajectories are depicted in Fig.l. 

Comparing with bare trajectories ai{t) in Fig.l one can see that the role 
and ordering of trajectories is changed as compared to bare ones when one 
goes from large t to the small t region. This property is very general and is 
not related to a particular choice of gik(t). 

It is also of interest to define the coupling of new Regge poles to the 
hadrons a, b and to probe in this way the quark and gluon contents of the 
poles. To this end we express the matrix Tik as 



ik 



where the diagonal matrix A is 



A 





V 







Oin^nO^k 





1 



(A5.13) 



j-Ol2{t) 










(A5.14) 



and find matrix elements Oik from the systems of equations 

aijOik + gi202k + gi^O-ik = 



(A5.15) 



'^The value of can be estimated from the model of /-dominance and experimental 

0.3 0.5. The yj"' coupling is 



(0) 



data on residues of the Pomeron and /-poles, 

(0) (0) 
512 513 ■ 



approximately ~ 0.5^^2^ while 523' 



36 



gi20ik + {ak - a2)02k + g2zOzk = (A5.16) 
g^Oik + g2i.02k + (ttfc - a'i)0^k = (A5.17) 
Equations (A5.17)-(A5.19) for A; = 1,2,3 and normalization condition 

|Ofcir + |Ofc2|' + |Ofc3|' = l (A5.18) 

define Oik up to a common phase. 

Physcially |Ofcjp = |Ojfcp gives a probability of finding original pole in 
the new pole k. 

Since original indices i refer to the glueball trajectory (i = 1), uu + dd 
/-trajectory {i = 2) and ss trajectory = 3), one can define in this way the 
percentage of the corresponding components in the new trajectory. 
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Table 1 

Effective mass eigenvalues fiQ{n,l) (in GeV for af — O.lSGeV^) obtained 

from Eq.(8), //q = ^o"adi(^^)^^'^- upper entry,and eigenvalues of reduced 

equation a(n)-lower entry. 



L\n 





1 


2 


3 





0.528 


0.803 


1.005 


1.174 




2.3381 


4.0879 


5.520 


6.786 


1 


0.693 
3.3613 


0.917 
4.8845 






2 


0.826 
4.2482 


1.020 
5.6297 







Table 2 

The eigenvalues (in GeV) of relativistic Hamiltonian for L = 



II 





1 


2 


3 


1 


5 




2.01 


2.99 


3.75 


4.37 


4.92 


5.41 



Table 3 

The eigenvalues (in GeV) of rotating string Hamiltonian (6) for L > 0. 



L 

n 


1 


2 


3 


4 


5 





2.65 


3.13 


3.53 


3.88 


4.206 


1 


3.645 


4.03 


4.366 


4.67 


4.95 


2 


4.40 


4.737 


5.04 


5.31 


5.56 


3 


5.02 


5.34 


5.62 


5.87 


6.10 


4 


5.58 


5.87 


6.13 


6.37 


6.59 


5 


6.09 


6.36 


6.60 


6.82 


7.03 
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Table 4 

Spin averaged glueball masses Mq/ ■yjo^ 



Quantum 

iiuiiil)ers 


This 

work 


Lattice data 


ref. [3] 


ref. [1] 


2 gluon 
states 


L = 0, = 


4.68 


4.66±0.14 


4.55±0.23 


L = 1, = 


6.0 


6.36±0.6 


6.1 ±0.38 


L = 0, = 1 


7.0 


6.68±0.6 


6.56±0.55 


L = 2, = 


7.0 


9.0 ±0.7(3++) 


7.7 ±0.4(3++) 


L — IjUr — 1 


8.0 




7.94 ±0.48 


3 gluon 
state 


K=0 


7.61 




8.19±0.48 



Table 5 

Masses of glueballs with L = 0, 1, 2 and n = 0, 1, a/ = O.lSGeV^ 



J'''' 


0-+ 


2+- 


0-+ 


2"+ 


n 





1 





1 





1 





1 


M(GeV) 


1.4 


2.4 


2.3 


3.3 


2.52 


3.55 


2.70 


3.7 





2++ 


0++ 


1++ 


3++ 


4++ 


3— 


2— 


1 — 


n 


























M(GeV) 


3.13; 3.11 


3.06 


3.07 


3.14 


3.16 


3.51 


3.23 


3.04 
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Table 6 

Comparison of predicted glueball masses with lattice data (for 

(7/ = 0.238G'e1/2 





M(GeV) 
This work 


Lattice data 


ref. [1] 


ref. [3] 


0++ 


1.58 


1.73±0.13 


1.74±0.05 


0++* 


2.71 


2.67±0.31 


3.14±0.10 


2++ 


2.59 


2.40±0.15 


2.47±0.08 


2++* 


3.73 


3.29±0.16 


3.21±0.35 


0-+ 


2.56 


2.59±0.17 


2.37±0.27 


0-+* 


3.77 


3.64±0.24 




2-+ 


3.03 


3.1±0.18 


3.37±0.31 


2-+* 


4.15 


3.89±0.23 




3++ 


3.58 


3.69±0.22 


4.3±0.34 


1 — 


3.49 


3.85±0.24 




2— 


3.71 


3.93±0.23 




3~ 


4.03 


4.13±0.29 
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Table 7 

Spin-orbit and tensor corrections to two-gluon glucball masses (the upper 
entries for n = and lower for n = l)af = O.lSGeV^ 







S=l 




L 


.=S=2 








0-+ 


2-+ 


0++ 


1++ 


2++ 


3++ 


4++ 




-0.197 


-0.0985 


- 0.1656 


- 0.138 


-0.083 





0.110 




-0.148 


0.074 


-0.128 


-0.107 


-0.064 





0.085 




-0.263 


0.0263 


-0.072 


-0.036 


0.015 


0.041 


0.020 




-0.198 


-0.02 


-0.056 


-0.028 


0.0116 


0.032 


0.0155 


^Kfrt 


-0.46 


+0.072 


-0.238 


-0.174 


-0.068 


0.041 


0.13 


-0.347 


+0.054 


-0.185 


-0.135 


-0.053 


0.032 


0.101 




0.082 


-0.041 


0.216 


0.18 


0.108 





-0.144 




0.05 


-0.025 


0.138 


0.115 


0.07 





-0.092 




-0.38 


+0.031 


-0.022 


+0.006 


0.04 


0.041 


-0.014 




-0.3 


+0.029 


-0.047 


-0.02 


0.017 


0.032 


0.009 


< Si2 > 


-2 


-1/5 


-2 


-1 


3/7 


8/7 


-4/7 


<LS> 


-2 


+1 


-6 


-5 


-3 
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Table 8 

Effect of inclusion of the Coulomb interactions on glueball masses and 
Regge paramerers, M(as, L = 0, 1, 2), cr/ = 0.18GeV^^ 



as 





0.2 


0.3 


0.39 


Mo{L = 0) 


2.11 


1.776 


1.587 


1.390 


Mo(L = 1) 


2.77 


2.56 


2.45 


2.36 


Mo(L = 2) 


3.30 


3.14 


3.05 


2.97 


«g(0) 


0.31 


0.298 


0.294 


0.290 




0.617 


1.06 


1.259 


1.44 
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Table 9 

Two-gluon glueball operators 



Dimension 


jPC 




^(m),(out) jj^ BPT 


4 


0++ 


tr{EiEi) 


tr{aiai) 


4 


2++ 


symrriik tr{EiEk) 


symrriik) tr{a,ai) 


4 


0-+ 


tr{EiBi) 


tr(aj(D X a)j) 


4 


1-+ 


tr(E X B) 


tr(a X (D X a)) 


4 


2-+ 


symm(^ik) tr{EiBk) 


symm(^ik) tr(ai(D x a)/,) 


4 


0++ 


tr{BiBi) 


tr(((D X a) ■ (D X a)) 


4 


2++ 


symnn^ik) tr{BiBk) 


symm{ik) tr{{D x a)j)(D x a)^) 


6 


3++ 


symm{kim) tr{D4^EkDiBm + D^BmDiEk) 


symm^iki) tr((D x a)i)Dfeai) 


6 


4++ 


symm^iikm) tr{DiEkDiEm) 


symm^ikim) tr{DiakDiam) 



Table 10 

Three-gluon glueball operators 



Dimension 


jPC 




\^(in),(out) jj^ BPT 


L 


BPT 


6 


1+- 




tr{{akai}(D x a),) 


1 


4 


6 


3+- 


symrukim tr{{Ek, Ei}B^) 


symmtr{{ak,ai}(D x a)„} 


1 


4 


6 


2+- 


symm^kn)enimtr{{Ek, Ei}B^) 


s|/mm(fc„)en/„ir({afc, ae}(D x a)^) 


1 


4 


6 


1— 


tr{EiEkEk) 


tr{aiakak) 





3 


6 


3— 


symm^kim) tr{EkEiEm) 


symm(kim) tr{akaiam) 





3 



45 



Table captions 

Table 1 

Effective mass eigenvalues Hoin, I) (in GeV for cr/ = O.lSGeV"^) obtained 
from Eq.(8), (Xq — upper entry,and eigenvalues of reduced 

equation a(n)-lower entry. 

Table 2 

The eigenvalues (in GeV) of relativistic Hamiltonian for L = 
Table 3 

The eigenvalues (in GeV) of rotating string Hamiltonian (6) for L > 0. 
Table 4 

Spin averaged glueball masses Mq / 
Table 5 

Masses of glueballs with L = 0, 1, 2 and n = 0, 1, a/ = O.lSGeV'^ 



Table 6 

Comparison of predicted glueball masses with lattice data (for af — 
0.238(761/2 

Table 7 

Spin-orbit and tensor corrections to two-gluon glueball masses (the upper 
entries for n = and lower ior n — l)af — O.lSGeV'^ 

Table 8 
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Effect of inclusion of the Coulomb interactions on glueball masses and Regge 
paramerers, M(as, L = 0, 1, 2), ct/ = O.lSGeV^ 

Table 9 

Two-gluon glueball operators 
Table 10 

Three-gluon glueball operators 
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0^ I 

1 2 3 4 5 

t 



Figure 1: Glueball, / and /' Regge trajectories as functions of t-channel 
energy squared t (in GeV"^). Dotted lines — bare trajectories, solid lines — 
trajectories with the coupling gik{t) taken into account in the form (A5.7), 
with parameters k — 1, = ^GeV'^. 
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